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ABSTRACT 

We study the spectroscopic properties of a selected sample of 26 events within Core Collapse 
Supernovae (CCSNe) family. Special attention is paid to the nebular oxygen forbidden line [O I] 
6300,6364A doublet. We analyze the line flux ratio "f63oo/f6364"> ^nd infer information about the 
optical depth evolution, densities, volume-filling factors in the oxygen emitting zones. The line lumi- 
nosity is measured for the sample events and its evolution is discussed on the basis of the bolometric 
light curve properties in type II and in type Ib-c SNe. The luminosities are then translated into oxy- 
gen abundances using two different methods. The resulting oxygen amounts are combined with the 
recovered ^^Ni masses and compared with theoretical models by means of the "[0/Fe] .v^. M,„s" 
diagram. Two distinguishable and continuous populations, corresponding to Ib-c and type II SNe, 
are found. The higher mass nature of the ejecta in type II objects is also imprinted on the [Ca II] 
729 1,7324 A over [O I] 6300,6364 A luminosity ratios. Our results may be used as input parameters 
for theoretical models studying the chemical enrichment of galaxies. 
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1. Introduction 

Recently, interest has increased the interpretation of Core Collapse Supernovae 
(CCSNe) data, both spectra and photometry (Richardson et al. 2006; Elmhamdi et 
al. 2006; Taubenberger et al. 2009; Maurer et al. 2010; Elmhamdi et al. 2011). 
In particular special attention has been devoted to the stripped-envelope events (i.e. 
type Ib-c hydrogen-deficient SNe). Studying enlarged samples of CCSNe objects, 
having good quality observations, can be a potential tool for assessing the similari- 
ties and the diversities within this SNe family, relating these facts to the physics and 
possibly to the nature of their progenitors. As an example, Elmhamdi et al. 2006 
have presented an investigation of the spectroscopic properties of a selected optical 
photospheric spectra of CCSNe, discussing how hydrogen manifests its presence 
within this class. The authors argued for a low mass thin hydrogen layer with very 
high ejection velocities above the helium shell to be the most likely scenario for 
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type lb SNe. Although the primary goal of the the cited work was highlighting the 
hydrogen traces in CCSNe, an important by-product result concerns the behaviour 
of oxygen, in particular the O I 7773 A Une. Based on the synthetic spectra fits, 
for this line it seems that at intermediate photospheric phases, type lb objects tend 
to have low optical depths, while some type Ic SNe, e.g. SN Ic 1987M, are found 
to display the deepest O I 7773 A profile. SNe of type lib & IT, at similar phases, 
are found to be the objects with the lowest 1 7773 A optical depth. At somewhat 
later epochs, transient type Ib/c objects display deep O 1 7773 A troughs. Mathe- 
son et al. (2002) arrived at similar conclusions, indicating the 1 7773 A line to be 
stronger in SNe Ic than in SNe lb. 

Interestingly, the deeper and stronger permitted oxygen line O I 7773 A in 
photospheric spectra of SNe Ic and Ib/c might imply that they are less diluted by 
the presence of a helium envelope. We expect indeed the oxygen lines to be more 
prominent for a "naked" C/0 progenitor core. Two further observational aspects 
tend to reinforce this belief: first, the forbidden lines, especially the [0 1] 6300,6364 
A doublet, seem to appear earlier following a SNe sequence "Ic — lb — lib — 11". The 
second indication comes from the fact that this nebular emission line has a velocity 
width decreasing following the SNe sequence above. 

In the present work we explore these points and other issues related to oxygen 
Unes within the CCSNe family, especially at late nebular phases. We analyze a 
sample of 26 CCSNe events, and quantify the [O 1] 6300,6364 A luminosities and 
discuss their evolution in time. The [Ca 11] 7291,7324 A over [O 1] 6300,6364 
A flux ratio is also shown and discussed as a possible progenitor mass indicator 
(Fransson & Chevalier 1989). Using the the computed luminosities we discuss 
two methods to translate these measurements into masses. The estimated oxygen 
masses are combined then with the measured nickel masses, determined from light 
curves, and compared with yields from three known theoretical CCSNe models by 
means of the oxygen-to-iron mass ratio. This is potentially important since our 
measurements can be directly correlated to the progenitor masses from theoretical 
models. 

Worth noting here is the importance of oxygen and iron estimates, from su- 
pernovae explosions, in the chemical enrichment and evolution of galaxies. In 
particular, the oxygen abundance is crucial in metal-poor stars, and is a key issue 
in modeling the early phases of the chemical galaxy evolution and as well in con- 
straining the age of globular clusters (VandenBerg et al. 2000; Melendez et al. 
2001). Oxygen is indeed considered as a major tracer of chemical evolution, since 
it is one of primary elements ejected by CCSNe (i.e. resulting from massive stars). 
For example the oxygen-to-iron ratio, when [0/Fe] > such as in the halo of our 
Galaxy, is an indicator of early chemical enrichment by massive stars. 

The paper is organized as follows. In Section 2, we describe the sample and 
we highlight the main spectroscopic characteristics and differences within the CC- 
SNe family. Some constructed quasi-bolometric light curves are also shown and 
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Fig. 1. Le ft panel : a sample of the late time CCSNe spectra. The most prominent lines are labeled. 
The corresponding observation date, since explosion time, are also reported. Right panel : the com- 
puted quasi-bolometric light curve of SN II 1987A (UBVRI bands); SN IIP 1999em (UBVRI-bands); 
SN lb 19901 (BVRI-bands); SN lib 1993J (UBVRI bands) and SN Ic I994I (BVRI bands). The Co 



to Fe decay slope is also shown for comparison. 



discussed. The [O I] 6300,6364 A line luminosity and the F6300/F6364 flux ratio 
are measured and presented in Section 3. In Section 4, we discuss two different 
methods for estimating oxygen mass in CCSNe. The methodology of ^^Ni mass 
estimation is given. We also evaluate the integrated flux ratio of the forbidden 
emission hnes [Ca II] 7291,7324 A and [O I] 6300,6364 A. We conclude with a 
summary and discussion of our findings in Section 5. 



2. The sample 

The selected sample consists of 26 CCSNe objects — 13 of them are type II, 
one of type lib and 12 of type Ib-c. Data are gathered mainly from the Uterature 
(i.e. pubUshed available data). Use is made of the "SUSPECT'Q and of the "CfA'l 
Supernova Spectra Archives. Some measurements are made on late spectra of SN 
1996aq, taken from the Padova-Asiago supernovae database. A summary of refer- 

' http://bruford.nhn.ou.edu/ ~ suspect/ indexl.html 

^ http : // cf a- w w w.harvard. edu/oir/Re search/supernova /SNarchi ve .html 
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ences and descriptive parameters of individual objects is given in Table 1. 

In what follows we adopt the standard reddening laws of Cardelli et al. (1989). 

An example of the analyzed late spectra is shown in Fig. 1 (left panel), together 
with identifications of the most prominent lines. We focus on the wavelength range 
5500-9200 A. At this epoch, when the events are in the radioactive tail phase, and 
except for the notable Ha emission in type II SNe, the optical spectra are domi- 
nated by emission lines of [O I] 6300,6364 A, [Ca II] 7291,7324 A and Ca II-IR 
triplet. The emission centered at X ~7800 A is usually attributed to O I 7774 A. 
These aspects will be highlighted when evaluating the [O I]/[Ca II] line intensity 
ratio of the sample objects. 

Worth recalling here that the [O I] 6300,6364 A line is found to be absent or 
very weak in SNe of type IIiH 

In Fig. 1 (right panel) the constructed "quasi-bolometric" light curves of SNIIP 
1999em, SNIIb 1993J, SNIb 19901 and SNIc 19941 are displayed. The pecuhar 
SNII 1987 A is also included for comparison. The available optical "U,B,V,R,I" 
broad-band photometric data have been integrated to recover the pseudo bolometric 
light curves. Literature references for the photometry, together with the adopted 
parameters are reported in Table 1 . We do not include the IR-contribution since it 
is available only for SN 1987A, and hence the derived integrated bolometric light 
curves represent a lower limit to the "real" bolometric light curves. The ^^Co to 
^^Fe decay slope, e— folding time of 111.3 days, corresponding to the full y— ray 
trapping is also shown. 

The difference in the early CCSNe bolometric light curves, i.e. the peak in 
type Ib-c SNe and the plateau in type II SNe, is mainly related to differences in 
presupernovae radii and structures. The plateau behaviour in type II is indicative 
of massive hydrogen envelopes, although its properties (duration and luminosity) 
predict also a dependence on radius, energy and the ejected amount of ^^Ni (Popov 
1993; Elmhamdi et al. 2003). The lack of such significant hydrogen in the outer 
layers of type Ib-c inhibits the presence of a plateau behaviour. 

The peak characteristics, luminosity and width, are sensitive to the ejecta mass, 
released energy and the ^^Ni mass (Arnett 1982). 

The clear faster decline at late phases for type Ib-c SNe is naturally attributed to 
the significant y-ray escape with decreased deposition as a result of low mass ejecta 
in this class of objects, while owing to the massive hydrogen mantle, type II SNe 
light curves indicate that radioactive decay of Co with the consequent trapping of 
y-rays is the main source of energy powering the light curves at late times. We note 

■'"n" stands for narrow. Type Iln are characterized by a narrow Ha emission and high bolometric 
light curve with a relatively flat evolution. In the prototype SN Iln 1988Z indeed there was no sign 
of the nebular forbidden lines [O I] 6300,6364 A and [Ca II] 7291,7324 A. Similar behaviour is 
observed in SN Iln 1994aj (Benetti et al. 1998). The well studied events SN Iln 1995N (Fransson 
et al. 2002; Pastorello et al. 2005) and SN UN 1995G (Chugai & Danziger 2003; Pastorello et al. 
2002) display the same characterizing peculiarity. One possibility for this is that their progenitors are 
not massive, believed to explode in a very dense CS environment. 
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that in type 11 objects the V-Ught curve follows the bolometric light curve fairly 
well at late epochs, which simplifies for example the derivation of the synthesized 
^^Ni mass. It is sufficient then the use of the tail absolute V-Ught curve of SN 
1987A , for which the ejected ^^Ni mass is accurately known from observations 
and detailed modeling, as a template for ^^Ni mass derivations in other 11 events 
(Elmhamdi et al. 2003b; Hamuy 2003). The stripped-envelope SNe case is more 
difficult. The V-band light curve does not parallel the bolometric one. It is hence 
necessary in type Ib-c a bolometric Ught curve modeling rather a simple use of the 
absolute visual bands. 
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Fig. 2. The temporal evolution of the [O I] 6300,6364 A luminosity for SNe of type II (upper panel) 
and for SNe Ib-c (lower panel). The Co to Fe decay slope is also reported. 
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3. The analysis 

In this section we measure and study the oxygen [O I] 6300,6364 A line lumi- 
nosity for the CCSNe sample. For this purpose the available spectra were corrected 
for redshift and reddening effects and calibrated with photometry when needed. 
The recovered integrated line fluxes, assuming a continuum level, are then con- 
verted to luminosities using adopted distances. 

3.1. Oxygen Luminosity and Mass 
3.1..1 The luminosity 

Figure 2 highlights the [O 1] 6300,6364 A line luminosity temporal evolution, start- 
ing at 200 days after explosion, for type 11 SNe (upper panel) and for type Ib-c SNe 
(lower panel). For comparison, the ^^Co to ^^Fe decay slope is also displayed (for 
an arbitrary ^^Ni mass; dashed line). 

The emission [O 1] 6300,6364 A light curves behave differently within the CC- 
SNe family. For type 11 events: the hght curves have a plateau-like maximum at day 
200, changing slowly until day ^ 340. At this epoch the hght curves are already 
on the exponential decline phase (Fig. 1 -right panel). The hne-luminosity dropped 
then sharply, with a rate of decrease similar to that of the radioactive decay of Co 
(i.e. e— folding time of 111.3 days). We note here that in the case of SN 1987A, it 
has been argued that dust condensation affects the hne luminosity evolution start- 
ing at day ~530, inducing a further increase in the rate of decrease (Danziger et 
al. 1991). At approximately the same time the [O 1] line profile showed a marked 
shift of the peaks towards blue wavelengths (Lucy et al. 1989). There are how- 
ever two exceptions decaying earher compared to the rest of the n SNe sample, 
namely SN 1970G and SN 1980K. Interestingly these two objects are respectively 
type IIP-L and IIL SNe. For type Ib-c events: starting at the age of 200 days, the 
light curves are already on a steep decline, faster to that of the ^^Co radioactive 
decay. The "Chi Square" fit to SN lb 1985F luminosity data for example indicates 
an e— folding time of about 70 days. The time at which the deviation to the decline 
occurs varies among the CCSNe family, being earlier in Ib-c, followed by IIL and 
then later on SNe IIP. Specifically the line luminosity trend is found to follow the 
bolometric hght curves in the time range of interest, namely after 200d in type Ib-c 
and after ~ 340d in type n. In SN 1987A for example, the [O I] luminosity relative 
to the bolometric one, i.e. {Lpj^/L^oi), shows an almost flat-topped behaviour at 
the time interval ~ 340 — 500 d. Other effects enter at later epochs, > 500d, such 
as dust condensation and the IR-catastrophe (Danziger et al. 1991; Fransson et al. 
1996; Menzies 1991). On the one hand this is a direct evidence that the dominant 
source of ionization and heating is y— rays from the radioactive decay of Co in 
the CCSNe variety, with the y— rays escaping with decreased deposition in type 
Ib-c events, owing to the low mass nature of their ejecta. We note also that in both 
panels the hne light curves span more than one dex in luminosity, which may be 
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related to the variation in the oxygen yields. Thus the importance of the oxygen 
mass estimates. 



[0 I] 6300,64 A at 390 d _ 




6200 6250 6300 6350 6400 
Wavelength (A) 



Fig. 3. The fit example to the [O I] 6300,6364 A profile of SN nP 1999em at day 390. Shown in 
thick line is the best fit. Dashed profiles refer to the decomposition of the best fit into two Gaussians 
for the [O I] 6300 A and [O I] 6364 A components. 



3.1..2 The flux ratio "Feaoo/FesM " 

Various investigators have dealt with the doublet flux ratio in the [O I] 6300,6364 
A line, especially for the well observed SN 1987A (Spyromilio & Pinto 1991; 
McCray 1996). It is well established that for a large and homologously expanding 
atmospheres such as in supemovae events, the Sobolev approximation holds and 
consequently simpUfies the radiation transfer and Une formation treatments (Castor 
1970; Jeffery & Branch 1990). 

Within this context, a hne intensity is given by: Flux <x p^^^ x A„/ x Nu, where 
Nu is the number density of atoms in the upper state, A„/ is the Einstein coefficient, 
and Pesc being the escape probability giving by Pgsc = [1 — exp{—i)]/i; i is the 
Sobolev optical depth. For a freely expanding atmosphere the density Nu goes as 

and hence the Sobolev optical depth should behave as i Nu x t = . 

For the Une of interest, i.e. the [O I] 6300,6364 A doublet, the two transitions 
have Agsoo = 3 x A^^^ , and T6300 = 3 x T6364 ; therefore the doublet flux ratio reads: 
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^6300 _ 1 -gxp(-T63oo) _ I - exp{-Xe3i.)o) 
^6364 i - exp{-i6364) I - exp{-ie3oo/^) 

A measure of the flux doublet ratio can be then translated into the optical depth 
in the 6300 A line in the SNe ejecta through solving Eq. 1 . In addition we expect 
an asymptotic value for the ratio to be 1 for the optically thick transitions and to be 
3 for the optically thin case. 

Recovering the ratio of the flux doublet is not always an easy task. Line blend- 
ing can prevent an adequate intensity measurement. One possibility is a simple 
use of the peak intensities in the components of the doublet assuming it similar to 
the ratio of the flux of the doublet components (Li & McCray 1992). Deblending 
yields to more accurate line ratio determination, especially in cases with a sever 
blend. We selected a variety of objects within our CCSNe sample that allow good 
estimates of the ratio. The studied events are: SNe IIP(1988A, 1999em, 2004et), 
SNe rb(1985F, 19901) and SNe Ic(1998bw, 2002ap). Data for SN 1987A are also 
included for comparison. An example is illustrated in Fig. 3 for the observed spec- 
trum of SN IIP 1999em at day 390. Our adopted methodology consists in fitting the 
[O 1] 6300,6364 A feature with a single function formed by 2 Gaussians, after esti- 
mating the continuum level and fixing the wavelength separation of the doublet, i.e. 
64A, and with both Gaussians having the same velocity width. The fit in Fig. 3 is a 
good one. In some cases it was not possible to obtain a particularly good overall fit 
using only two Gaussians since one is left with either a residual in the blue wing of 
the 6300 A component for narrower lines, or Uttle evidence of the the 6364 A com- 
ponent for broader lines. Elmhamdi et al. 2004 have used similar methodology and 
obtained a satisfactory fit in the specttum of SN lb 19901 at 254d inttoducing a third 
Gaussian component assuming there is a contribution from Fe IT 6239 A emission, 
and allowing its velocity width to be a free parameter (their Fig. 7). According 
to our sample analysis the two approaches, using the deblending technique and a 
simple use of the peak fluxes, when possible, agree within 20 % . Figure 4 reports 
the results of our investigation. The upper panel displays the temporal evolution of 
the ratio measurements, while in the lower panel we plot the corresponding optical 
depths in the primary component, computed by solving Eq. 1. The ratio values 
demonstrate a temporal trend from an optically thick phase to an optically thin one. 
The optical depth is found to decrease from values as high as 7 around day 170 
falling to lower values at later phases (< 0.5 around 500d). The continuous line is 
the resulted LSQ fit to the data, using a power-law function (i.e. <=< t^"). The best 
fit gives an index of n= 2.089 (±0.153). This is consistent with the index expected 
from the simple expansion assumption discussed above (i.e. T<x NuXt = t~^)- 

Furthermore, we may introduce a time dubbed ttrms > as the time at which the 
line makes the transition to the phase characterized by T6300 = 1 . In fact using 
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Fig. 4. Upper panel: the temporal evolution of the flux ratio ^6300 7^6364 • The overploted curve 
corresponds to the LSQ-fit by Eq. 2, indicating a value of tnans — 426d (see text for details). 
Lower panel : the recovered optical depths in the [O I] 6300A line. The power law LSQ-fit is also 
displayed (^J"" with n = 2.089 ±0.153; see text). 

T o< and introducing ttrans > Eq. 1 can be written as: 

^6300 ^ ^ - exp{- {ttrans /t f) 
^6364 1 - exp{- {ttrans /V^t)^) 

Provided the ratio measurements, we fit the data with Eq.2. The LSQ fit is shown 
by the continuous line in the upper panel of Fig. 4. The best fit indicates a time of 
ttrans = 426.06 it 13.24 d. The time at which T6364 = 1 , is 245.96 ±7.64 d. 

This above described analysis demonstrates the consistency within the CCSNe 
family on how the components intensity ratio and the optical depths develop in 
time. 
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3.1..3 The mass 

In this section we adopt two different methods for estimating the oxygen mass in 
SNe of type II and SNe of type Ib-c. 

• Type 11 SNe: 

We use the recovered [O I] doublet luminosities to determine the abundance of 
the oxygen produced through the SN explosion. 

At the epoch of about 1 year, the luminosity of the [O I] doublet is powered 
by the y-ray deposition and by ultraviolet emission arising from the deposition of 
y-rays in oxygen material. The [O I] doublet luminosity is related directly to the 
oxygen mass (Elmhamdi et al. 2003a), and at a given time one may write: 

L,([OI])=TixL,(56Co)x^ (3) 

where Mq is the mass of oxygen, Mgx is the "excited" mass in which the bulk 
of the radioactive energy is deposited, and r\ is the efficiency of transformation 
of the energy deposited in oxygen into the [O I] doublet radiation. The L((^^Co) 
refers to the radioactive decay energy input, given by = LqX {MiiilMQ)Q~'l'^° , 
with Lo ~ 1.32x 10^^ ergs s~^ , the initial luminosity corresponding to 1 Mq and 
Tco = 111.3 d. 

Assuming then that all type II events have similar r\ and Mgx at similar phases, 
we derive rough estimates of the oxygen mass for the events of the sample adopting 
[O I] 6300,6364A Hght curve of SN 1987A as template. As discussed previously, 
the [O I] luminosity is found to have similar decay rates in SNe II at late phases, 
following a behaviour similar to the bolometric light curves. This fact gives some 
confidence for the use of SN 1987A [O I] luminosity as a template for the recovery 
of the amounts of ejected oxygen. The oxygen mass in SN 1987A is estimated to 
be in the range \.5-l Mq (Fransson et al. 1996; Chugai 1994). 

In Table 1, column 7, we report the amounts of the ejected oxygen mass derived 
in this maimer. The variation range, for each event, reflects the combination of the 
variation in the oxygen mass for SN 1987A (i.e. 1.5 — 2 Mq) together with the 
uncertainties from the fit procedure. 

• Type Ib-c SNe: 

The emission lines which are formed at densities above their critical ones, given 
the optical depths are not large, have the luminosity directly proportional to the 
mass of the emitting ion. Such conditions hold for the [O I] 6300,6364 A doublet 
line emission in type Ib-c SNe at nebular phases. On the one hand, we found 
moderate optical depths at late phases (Fig 4, lower panel). On the other hand, 
the condition of the high density limit above the critical density for the [O I] line 
(~ 7 X 10^ cm~^) is found to be clearly fulfilled in the ejecta of type Ib-c SNe 
(Leibundgut et al. 1991; Elmhamdi et al. 2004). A possible direct method can be 
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used to check the high density limit characteristic. Indeed, the density is directly 
related to the relative strengths of the [O 1] doublet components (Leibundgut et al. 
1991; SpyromiUo & Pinto. 1991). In both cited works, the variation of the line 
doublet ratio as function of the density at a given time is computed, and is found 
to be insensitive to the adopted temperature, especially for the late epochs. From 
our computations we take a ratio value of 1.58 at day 250 or 1.8 at day 300 as 
representatives (i.e. from the best fit in Fig.4-upper panel). According to Fig. 6 of 
Leibundgut et al. (1991), the uncertainty in the temperature leads to the following 
density range 2 x 10^ cm^^ < A'^^, < 4 x 10^ cm^^* . 

In the high density limit, i.e. above the critical density, the mass of ejected 
oxygen can be recovered using the [O I] 6300,6364 A flux. Uomoto (1986) has 
shown that the oxygen mass, in Mq , is given by: 



Mox = 10^ X Z)2 X F{[0 I]) X exp{2.28/T4) (4) 

where D is the distance to the supernova (in Mpc), F is the reddening-free [O I] 
integrated flux (in ergs s~^ cm ~^) and is the temperature of the oxygen-emitting 
gas (in 10^ K). 

Worth noting that because of the variation of F([0 I]) and 74, Eq. 4 implies 
time-dependence. Schlegel & Kirshner (1989), when estimating the ejected oxy- 
gen amounts, have adopted a constant temperature 74=0.4 at the nebular phase of 
SNe lb 1984L and 1985F. The assumption of a constant temperature at different late 
phases provides different oxygen masses as one may expect, since earlier nebular 
epochs are hotter compared to latter ones. 

Alternatively, the temperature at a given time can be constrained based on the 
[O I] 5577 A to [O I] 6300,6364 A flux ratio. Assuming that the O I hues are 
formed mainly by coUisional excitation, the ratio is given by the following expres- 
sion (Houck & Fransson 1996): 

= 0.03P6300 X [1 + 1.44r,-^''^^(:^)] X exp{25.83/Ts) (5) 



where Tj, is the temperature of the oxygen-emitting gas (in 10^ K), Pesoo is the 
[O I] 6300 A Sobolev escape probability (~ 1) and Ne is the electron density. 

A problem with this method is the very weak observed [O I] 5577 A feature in 
the late time spectra analyzed here. It is indeed an indication of low temperatures 
in the oxygen-emitting zone. However one may estimate an upper limit on the [O 
I] 5577 A flux by integrating over the same velocity interval of the [O I] 6300,6364 
A profile. It is found that the temperatures at late phases of interest, > 250d, tend 
to have values in the range 3400 — 4200 K. 

Results from this described method, through equation 4, are listed in column 
7 in Table 1. For SNe 19941, 1993J, 1998bw and 2002ap however, the amounts 
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Fig. 5. The temporal evolution of the [Ca II] 7291,7324 A over [O I] 6300,6364 A integrated flux 
ratio. Note the complete separation of type Ib-c from type II SNe. 

reported in the table come from the most recent spectroscopic and photometric 
modeUng (see column 8 for references). 

3.1..4 The [Ca II]/[0 I] intensity ratio 

After considering the reddening correction, we have evaluated the integrated flux 
ratio of the forbidden emission lines [Ca 11] 7291,7324 A and [O 1] 6300,6364 A, 
for the CCSNe sample objects having wavelength coverage in these two nebular 
lines. Results are displayed in Fig. 5. Two separate populations are clearly distin- 
guishable. The ratio is found to remain below unity in type Ib-c, with the exception 
of SN 19901b around day 315 (~ 1.22). SNe Eb 1993J, Ic-hypemova 1998bw and 
ILL 1 980K belong to this category as well. A mean value is measured to be ~ 0.5 1 . 
SNe of type n instead concentrate on the top of the figure, with a mean value of 
~ 3.17. Additionally, the two classes appears to have a flat evolution behaviour of 
the [Ca ll]/[0 1] ratio. 

In a detailed analysis, Fransson & Chevalier (1987, 1989) have modeled late 
emission spectral lines for different supernova progenitors. It is found that be- 
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cause of the composition and density structures one can use the relative emission 
line strengths as a progenitor indicator tool. In particular the forbidden emission 
line ratio [Ca n]/[0 I] presents weak dependence on density and temperature of 
the emitting zones, and is expected to display an almost constant evolution at late 
epochs. The ratio is found to be very sensitive to the core structure and mass. 
Furthermore it seems that the ratio tends to increase with decreasing progenitor 
mass. The distribution of our measurements in two groups tend to indicate differ- 
ent progenitor properties, with lower progenitors for type Ib-c, lib and IIL classes 
compared to normal type IIP events. It is worth to note here however that in type 
Ib-c SNe there is no hydrogen rich Ca II emitting zone as is the case for type n 
objects (de Kool et al. 1998). 

3.1..5 The 56 Ni mass 

As discussed in Sect. 2, the exponential behaviour of the late V-band absolute 
light curves of type 11 SNe is found to be in accordance with the radioactive de- 
cay 5^Co— 5'^Fe (My oc exp{—t/l\\3 d)), which provides a potential method to 
recover the amount of the ejected ^^Ni with the use of the My photometry of SN 
1987A as a template between 120-400 d (Elmhamdi et al. 2004; Hamuy 2003). 
An amount of 0.075 Mq is adopted for the ^''Ni mass of SN 1987 A (Catchpole 
et al. 1988, Suntzeff & Bouchet 1991). Results from the described methodology 
are reported in Tab. 1 (column 6). There is a significant scatter of the ejected ^^Ni 
masses, with an average value of « 0.053 Mq . 

For type Ib-c events, the broad band light curves do not trace necessary the 
bolometrics. The amounts summarized in Tab. 1 (column 6) are indeed results from 
bolometric Ught curves and/or spectra modeUng (see column 8 for the correspond- 
ing references). An average value of 0.18 Mq is computed. 

4. Discussion and conclusions 

We have selected a sample of 26 events within the CCSNe family. Our main 
goal was to investigate how the oxygen manifests its presence at late phases for 
each SN class, especially the emission [O I] 6300,6364 A doublet. 

Based on investigating early spectra of CCSNe, the permitted oxygen Une 
O I 7773 A seems to get weaker following the SNe sequence "Ic — lb— lib— 11" 
(Elmhamdi et al. 2006; Matheson et al. 2001). For stripped envelope objects, be- 
ing less diluted by the presence of an helium envelope, one may expect the oxygen 
lines to be more prominent. An observational fact that clearly support this belief 
is the earUer appearance of the [O I] 6300,6364 A following the above order, i.e. 
"Ic — lb — lib — 11". The following examples highlight this fact. The oxygen line 
emerges at an age of 1-2 months in type Ic SN 1987M (Filippenko 1997). SN Ic 
19941 displayed evidence for the line at an age of 50 days, although some hints may 
even be seen in the ~36 days spectrum (Clocchiatti et al. 1996b). SN Ic 1997B 
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shows clear evidence for [0 1] 6300,6364 A, [Ca 11] 7291,7324 A and also 1 7773 
A on the 2 months spectrum (Gomez & Lopez 2002). In SN 1998bw, Ic-hypemova 
event, the nebular features were already recognizable in the 43 days spectrum (Patat 
et al. 2001). While in SN lb 19901 it was hinted at the 70 days spectrum (Elmhamdi 
et al. 2004). In other type lb SNe it appears earlier than in SN 19901. In SN lib 
1993J, a transition object, the hne was visible in the 62 days spectrum (Barbon et 
al. 1995). SN 1996cb, another well observed lib event, showed evidence of the [O 
1] 6300,6364 A line around day 80 (Qiu et al. 1999). In SNe 11, however, the line 
appears later: around day 150 in SN 1987A (Catchpole et al. 1988) and after day 
138 in SN 1992H (Clocchiatti et al. 1996a). In SN 11 1999em it is suggested at a 
somewhat earUer phase compared to SNe 1987A and 1992H, namely at day 114. 

Furthermore, the width velocities in the [O I] 6300,6364 A and in the [Ca 11] 
7291,7324 A nebular lines are mainly found to be higher in SNe Ic than in SNe 
lb (Schlegel & Kirshner 1989; Matheson et al. 2001). In the cited papers, the 
FWHM values were evaluated fitting a single profile, namely Gaussian, to the total 
observed profiles. We mention here that given the line blending in these two Unes 
(doublets), results should be interpreted with some caution. One should indeed 
resolve the lines and compare width velocities of single components rather than 
using a single Gaussian fit to the whole feature. It is not simple however in some 
cases to get a satisfactory fit, especially when the blend is sever. According to our 
analysis, we faced more difficulties in type lb and Ic than in type 11. This fact is 
indeed in itself indicative of higher velocities in stripped envelope objects, while 
in SNe II, owing to lower expansion velocities, the observed [O I] 6300,6364 A 
profile is less complex and the two components are clearly visible. This can be 
easily understood if considering the formation of the line doublets in cases where 
the element expansion velocity does not exceed the the two components separation 
velocity, i.e. 3000 km s~^ . We have checked the width velocities by means of the 
FWHM of the single [O I] component at 6300 A of some of our sample events. 
At a phase of '^300d, SNe 11 1987A and 1999em are found to have respectively 
FWHM {[O I] 6300 ) ~ 2750 and 2400 km s^^ , while for lb SNe a typical range 
variation of ~ 4000-5000 km s"! is deduced (e.g. SNe lb 1985F, 1996N). SN 
nb 1993 J, at similar age, has a velocity of ~ 4600 km s~^ , with a complex [O I] 
6300,6364 A profile due to the presence of Ha at the red wing of the line (Patat 
et al. 1995). Higher velocities, > 6000 km s^', are computed from Ic spectra 
(e.g. SNe 1987M, 19941 and 1998bw). The [O I] 6300,6364 A nebular hne at late 
epochs is representative of the expansion velocity. For a given explosion energy, a 
greater ejecta mass allows for lower velocities. The found lower velocities in type 
n SNe can be attributed then to a large mass with respect to type Ib-c objects. 

We have investigated in details the properties of the [O 1] 6300,6364 A doublet 
profile. The measured line luminosities are found to trace well the bolometric hght 
curves, behaving differently in type Ib-c than in type n events. In SNe n, the hne 
hght curves show a plateau-hke evolution until ~ 340 days, following later-on by 
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an exponential decline similar to the radioactive decay of Co (e—folding~ 111.3 
days). The linear type n SNe display an early decline in their [O 1] light curves, 
type Ib-c light curves display already a decUne by an age of 200 days, steeper than 
the ^^Co radioactive decay rate. The rays from the radioactive decay of ^^Co 
are hence the dominant source of ionization and heating in CCSNe variety. Further- 
more, owing to the low mass nature of their ejecta, type Ib-c SNe are characterized 
by a decreasing deposition of the y— rays escape. Similar conclusions about the 
progenitor diversity are argued from the analysis of the integrated flux ratio of the 
forbidden emission lines [Ca 11] 7291,7324 A and [O 1] 6300,6364 A. The ratio is 
found to concentrate in two distinguishable locations, namely around a value of 3 
for type n SNe (mean ~3.17), and below unity in Ib-c objects (mean ~0.51). This 
forbidden lines ratio is potentially sensitive to the progenitor mass star (Fransson 
& Chevalier 1987; 1989), suggesting higher masses in SNe of type 11. 

The way CCSNe family members transit from optically thick to optically thin 
phases is imprinted onto the profile of the [O 1] 6300,6364 A doublet. We empha- 
size this point making use of the components flux ratio F63oo/^6364^ after decon- 
volving the two components of the observed line profiles. Based on a simple de- 
scription of the physics behind the line formation, our results indicate a consistency 
within the CCSNe on how the components intensity ratio and consequently the op- 
tical depths develop in time. Additionally, the ratio F63oo/^6364 can be used to 
estimate the average density of the oxygen-emitting zone (Leibundgut et al. 1991; 
Spyromilio & Pinto 1991). Adopting representative values of the ratio of 1.58 at 
day 250 or 1.8 at day 300, inferred from our best fit in Fig.4-upper panel, and ac- 
cording to Fig. 6 of Leibundgut et al. (1991), the uncertainty in the temperature 
leads to the following density range 2x10^ cm" <Ne<4x 10'' cm . Worth 
noticing here that based on our approach one may get further informations about 
the physical conditions in the oxygen emitting region. In fact, for an uniform den- 
sity distribution, one may use the above derived density range together with the 
element volume, based on the [O I] 6300 A FWHM velocity, to get a rough es- 
timate of the ejected oxygen mass. This approach is found to provide too large 
amounts compared to results reported in Tab. 1 (column 7). For SN IIP 1999em 
for example, at 300 days, using the previous derived velocity of 2400 km s~^ , we 
find an ejected oxygen mass as high as ~40 Mq , indicating a volume-filling factor 
as low as lO^^. Similar filhng factor is found for SN 11 1987A and SN lb 19901. 
These results suggest that the oxygen material in CCSNe generally fills its volume 
in a clumpy way rather than homogeneously. 

We describe and adopt two methods for the oxygen mass estimate in CCSNe. 
We provide also estimates of the ejected ^^Ni masses. It should however be noted 
that various parameters affect the recovered values, especially as long as the oxygen 
mass is concerned, contributing hence to their uncertainties. Indeed, the luminosity 
and temperature measurements are sensitive to the nature and quality of the spectra 
(i.e. severe line blending; multi-component profiles; presence of [O I] 5577 A). 
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According to our test effects, we estimate the uncertainty of the derived oxygen- 
mass amounts reported in Tab. 1 to vary within 15 — 30%. Furthermore, for cases 
with clear evidence of O I 7774 A line, the presence of ionized oxygen is argued 
(Mazzali et al. 2010), requiring higher mass of oxygen than the one needed to pro- 
duce [O I] 6300,6364 A alone. Generally, Nickel masses are found to suffer less 
uncertainties compared to oxygen. The errors in type II SNe, being more homo- 
geneous than SNe Ib-c, do not exceed 20% (mainly due to distance and reddening 
estimates). In SNe Ib-c, we estimate a variation of 10 — 30% ." 

These estimates , i.e. oxygen and iron masses, are of importance as they can be 
indicative of the core mass and related to the progenitor star nature. For the purpose 
of a better investigation, we compute the [0/Fe£| yields ratio for each event. A 
solar value of logio{0/Fe)Q = 0.82 dex (Anders & Grevesse 1989) is adopted. 
The recovered amounts are reported in Fig. 6 (horizontal lines) as function of the 
initial mass according to the most reliable core collapse SN models in literature , 
namely Woosley & Weaver 1995 (short dot-dashed line), Thielemann, Nomoto & 
Hashimoto 1996 (long dashed line) and Nomoto et al. 1997 (long dot-dashed line). 
A typical value for type la SNe is as well reported (Nomoto et al. 1984), which 
shows the nature of type la SNe being iron producers. 

Fig. 6 defines two possible concentration regions, of type II objects from one 
side and type Ib-c from the other side, and suggestive of a continuum in the [0/Fe] 
values. Indeed, an important issue that can be read out immediately from Fig. 6 is 
that the reported results for type Ib-c are found to be located at the bottom of type 
II SNe. These zones may provide constraints on the progenitor masses in CCSNe 
family. Within the core collapse models of type Ib-c SNe, two scenarios are argued: 
first a single high mass star {M^s > 35 M0) exploding as Wolf-Rayet star after 
an episode of strong stellar wind, and second a less massive star (M,„v 13 — 18 
Mq ) in a binary system. Although the progenitor nature of type Ib-c SNe is an open 
and debated issue, the position of type Ib-c events in the "[0/Fe] .vs. M,„s" plot, 
according to Fig. 6, might be taken as an observational support of the intermediate 
massive stars in binary systems, stripped of their envelope through binary transfer, 
as the favoured progenitors in this class of objects. 

Although the uncertainties in the stellar evolution models of massive stars and 
in the determination of the oxygen and iron yields. Fig. 6 provides a methodology 
to elucidate the progenitor nature of core collapse SNe and interesting comparisons 
can be drawn . However, extended samples and more reliable determinations, es- 
pecially for oxygen abundances (late spectra modeling for example), are clearly 
needed to populate the "[O/Fe] .vs. M,„v" diagram and hence have a deeper view 
on the oxygen to iron ratio and how it changes as function of the type of SN. The 
present quantitative analysis can also provide insights and input data for the evolu- 
tion of galaxies and the chemical enrichment models. 

The main goal of the paper was to highlight some late-phases spectra related 
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properties of CCSNe, providing possible physical interpretations. As a first step 
we concentrate our study on 26 selected events. For future investigation we aim 
to enlarge the studied sample including more recent and up-dated observed CC- 
SNe objects, especially those with spectra of improved time coverage and spectral 
resolution (mostly last decade observations; Elmhamdi et al. In preparation). 
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Fig. 6. The [0/Fe] yields ratio as a function of the initial mass according to the theoretical mod- 
els of Woosley & Weaver 1995 (short dot-dashed line), Thielemann, Nomoto & Hashimoto 1996 
(long dashed line) and Nomoto et al. 1997 (long dot-dashed line). For each event the corresponding 
estimated [0/Fe] value is plotted horizontally (SNe 11 with horizontal dashed Unes, while horizon- 
tal dotted lines correspond to SNe Ib-c). Typical value for thermonuclear SNe is also reported for 
comparison (horizontal continuous line, Nomoto et al. 1984). 
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Table 1 : Parameters data of the CCSNe sample 



SN 


SN 


Parent 


Distance 




Mm 




References 


name 


Type 


galaxy 


(Mpc) 




Mq 


Mq 




1987A 


n 


LMC 


0.05 


0.6 


0.075 


1.5-2 


1,2 


1988A 


np 


NGC 4579 


22.95 


0.136 


0.088 


0.7-1.38 


7,8 


1988H 


IIP 


NGC 5878 


28.47 


0.47 


0.033 


0.79-1.1 


8 


1990E 


IIP 


NGC 1035 


16.18 


1.2 


0.043 


0.94-1.3 


10, 13 


1991G 


np 


NGC 4088 


14.11 


0.065 


0.021 


0.36-0.49 


9 


1992H 


HP 


NGC 5377 


29.07 


0.054 


0.123 


0.98-1.5 


14 


1995ad 


IIP 


NGC 2139 


23.52 


0.112 








17 


1997D 


IIP 


NGC 1536 


16.84 


0.07 


0.0065 





11, 12, 15 


1999em 


np 


NGC 1637 


8.8 


0.31 


0.027 


0.3-0.54 


3,4,5 


2002hli 


HP 


NGC 6946 


5.9 


3.3 & 1.7 


0.07 


0.78-1.05 


18 


2004et 


IIP 


NGC6946 


5.6 


1.27 


0.06 


0.75-1.8 


19 


1970G 


IIL-P 


NGC 5457 


7.2 


0.4 


0.051 


0.88-1.18 


6, 16 


1980K 


nL 


NGC6946 


5.6 


1.5 


0.046 


0.57-0.76 


35, 36 


19901 


lb 


NGC4650 


39.3 


0.4 


0.11 


0.7-1.35 


20 


1996aq 


Ib/c 


NGC5584 


24.32 


0.124 




0.4-1 


17 


19831 


Ic 


NGC4051 


22.08 


0.043 


0.15 


0.43 


21 


1983N 


lb 


NGC5236 


4.46 


0.51 


0.15 


0.43 


21,22 


1990B 


Ib/c 


NGC4568 


15.27 


2.64 


0.14 


0.4 


23, 24, 25 


1984L 


lb 


NGC991 


23.44 


0.32 


0.37 


0.58 


25, 26 


1985F 


lb 


NGC4618 


7.3 


0.0 


0.1 


0.3 


26 


1987M 


Ic 


NGC 2715 


22.7 


1.3 


0.26 


0.4 


27, 28 


19941 


Ic 


NGC 5194 


8.32 


0.93 


0.07 


0.22 


29 


1993J 


lib 


M81 


3.64 


0.6 


0.1-0.14 


0.5 


30, 31, 32 


1996N 


lb 


NGC 1398 


22 


0.0 




0.1-0.3 


30 


1998bw 


Ic 


ESQ 184-G82 


35.16 


0.2 


0.4-0.5 


5-6 


33 


2()()2ap 


Ic 


M74 


7.9 


0.24 


0.09 


0.6 


34 
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NOTE: the net effect from the uncertainty sources, see text, indicate a variation in the calculated oxygen masses 
within 15 - 30% . The estimated iron mass is found to be less biased, especially for type 11 events. 



